X-ray diffraction measurements of the c-axis Debye- Waller factors of YBa 2 Cu 3 07 and 

HgBa 2 CaCu 2 6 
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We report the first application of x-rays to the measurement of the temperature dependent Bragg 
peak intensities to obtain Debye-Waller factors on high-temperature superconductors. Intensities 
of (0,0,/) peaks of YBa2Cu3C>7 and HgBa2CaCu2C>6 thin films are measured to obtain the c-axis 
Debye-Waller factors. While lattice constant and some Debye-Waller factor measurements on high 
T c superconductors show anomalies at the transition temperature, our measurements by x-ray 
diffraction show a smooth transition of the c-axis Debye-Waller factors through T c . This suggests 
that the dynamic displacements of the heavy elements along the c-axis direction in these compounds 
do not have anomalies at T c . This method in combination with measurements by other techniques 
will give more details concerning dynamics of the lattice. 

PACS numbers: PACS numbers: 74.25.-q, 61.10.Eq, 74.25.Kc 
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The majority of efforts to understand the mechanism 
of the high T c superconductivity (HTSC) have been in 
terms of pure electronic effects, primarily because of the 
conjecture that electron-phonon coupling alone cannot 
give such a high T c . However, there have been persistent 
efforts to study the role of the lattice in HTSC 1 : 2 : 3 : 4 ! 5 ! 6 ! 7 . 
In addition to the fact that the structural information 
provides important information about the materials, the 
lattice can still play an important role as it affects both 
the hopping (t) and magnetic exchange (J) energies of 
these materials, the two most important parameters in 
the physics of these correlated systems. Indeed, there 
have been theoretical efforts to study the role of lattice 
vibration in the cupratesi*^. Therefore, experimental in- 
vestigation of the temperature dependent lattice vibra- 
tion in these materials, especially the critical behavior 
around T c , is essential and can provide vital information 
concerning the low energy degrees of freedom. 

Experimental evidence that shows clear correlation be- 
tween the critical behaviors in the electronic and crys- 
tal structures are anomalies in lattice constants at T c . 
These have been observed in various HTSC materials by 
x-raji&4 and more accurately by capacitance dilatomctry 
method a 5 ! 6 : 7 . They have shown that the lattice constant 
generally decreases at a much faster rate at T c than at 
other temperatures as the temperature decreases. Sub- 
sequently, the thermal expansion coefficient a shows an 
anomaly that is of similar shape to the heat capacity 
anomalies observed at T c . This enhanced decrease in the 
lattice constant has been related to the specific heat ca- 
pacity through the thermodynamic Ehrenfest relations 7 . 
Yet the microscopic reason behind such a lattice anomaly 
is not well understood. One natural question is if it is re- 



lated to the dynamic properties of phonons. It will there- 
fore be interesting to measure the Debye-Waller factors 
(DWF) of HTSC materials. 

In complex materials, the DWFs are usually measured 
by extended x-ray absorption fine structure (EXAFS). 
In the HTSCs there have been efforts to measure the 
temperature dependent mean-square displacement a 2 (T) 
(which is inversely related to the DWF) by EXAFS 8 ' 9 ' 10 , 
with essentially only Hg based samples showing anoma- 
lies at TfiH. EXAFS measures local coordinates at a 
certain atom and therefore has an important advantage 
of being element specific. It, however, measures only 
relative atomic motions (for example, O relative to Cu 
or optical phonons in the Cu-0 plane) and is insensi- 
tive to collective motions of the atoms, for instance, the 
low frequency acoustic phonons. To see the effect in the 
phonon population, the necessary condition is to have 
ample phonons in the first place. This may be difficult for 
optical phonons as their energy scale is usually very high 
compared to the experimental temperature scales (which 
are set by T c ). This may be related to the fact that 
EXAFS results largely show no anomalies. In contrast, 
even though they have their own disadvantages which 
will be discussed later, diffraction techniques in principle 
measure the integrated phonon effects^. Here we report 
novel application of x-ray diffraction that can in principle 
measure DWFs. This is to our knowledge the first ap- 
plication of this method by x-rays to complex materials 
like HTSCs. We have measured the temperature depen- 
dence of the Bragg peak intensities, a quantity that can 
be related to the dynamic lattice motion. The results are 
discussed along with various aspects of x-ray diffraction 
and possible future experiments. 



The experiments were performed at beamline 2-1 of 
the Stanford Synchrotron Radiation Laboratory (SSRL) 
which is equipped with a 2-axis Huber diffractometer. 
The diffractometer has an open cycle He cryostat. hv — 
8.8 keV x-rays, just below the Cu K edge, were used to 
reduce the background from the fluorescent light. For the 
reasons explained below, intensities of the (0, 0, 1) peaks 
of YBa 2 Cu 3 7 (YBCO) and HgBa 2 CaCu 2 6 (Hgl212) 
thin film samples grown on SrTiC>3 substrates were taken. 
These materials were grown according to the procedures 
described elsewhereiSJ^. The YBCO film has T c of 
90. 5K with the AT C less than IK. The Hgl212 sample 
has the AT C less than IK locally but shows distribution 
of T c 's from 123K at the center to 118K at the edge. 

In order to obtain highly reliable intensities, it is es- 
sential to reduce the errors from the beam and sample in- 
stabilities and detector non-uniformity. The beam drifts 
were monitored in the Io section and corrected for before 
each measurement. The biggest effect came from the mo- 
tion of the cryo-head during the temperature cycle. To 
reduce such effects, a special manipulator has been de- 
signed as shown in Figure 1. The sample holder is in weak 
thermal contact with the cryo-head and has its own tem- 
perature sensor and heater. The thermal contact is weak 
enough so that the temperature of the sample holder can 
be varied between 10 K and 200 K with the temperature 
of the cryo-head maintained at 10 K. This removed most 
of the sample motion and greatly improved the repro- 
ducibility of the data. 

The special manipulator did not completely remove the 
sample motion with respect to the beam during tempera- 
ture cycling. If the sample was completely homogeneous 
then the relative motion of the sample will have no ef- 
fect on the measurement. In our experience, the only 
high T c samples which are homogeneous to a highly col- 
limated and monochromatic synchrotron beam are fine 
grain powders. However, in a powder sample the high 
index peaks, the peaks more sensitive to the DWF (as 
DWF ~ I 2 ), were too weak to measure with sufficient 
statistical accuracy in a reasonable amount of time. Even 
though the intensities of the high index Bragg peaks were 
acceptable in a single crystal we discovered, after looking 
at high T c single crystals from several different sources, 
that these samples were really a bundle of several low an- 
gle grain boundary crystallites. Even a small motion of 
such a sample with respect to the incident beam changed 
the fraction of the diffracting domain in the beam, which 
resulted in a significant change in the diffracted beam in- 
tensity. Thin films of sufficiently high quality and lateral 
homogeneity for synchrotron measurements were avail- 
able but they were all (0, 0, 1) orientationM. Hence, as 
a compromise between high intensity and sufficiently ho- 
mogeneous samples, we performed our measurements on 
the (0, 0, 1) peaks of thin film samples. 

In addition to the motion of the cryostat upon the 
temperature change, other factors such as detector non- 
uniformity contributed to the irreproducibility of the 
data. Since the effect we were looking for was very small, 




FIG. 1: (a) The schematic of the cryostat designed for the ex- 
periments. The sample holder is in weak thermal contact with 
the cryo-head with its own temperature sensor and heater. 
Due to the weak thermal contact, the temperature of the sam- 
ple holder can be varied between 10 K and 200 K while the 
temperature of the cryo-head fixed at 10 K. (b) Illustration 
of the experimental configuration. The detector is fixed at 
the 29 position while 9 is scanned. The detector slit is wide 
enough to accept the whole Bragg peak as illustrated unlike 
the conventional setup. The 29 position is adjusted at each 
temperature so that the same part of the detector is used. 



the accuracy/reproducibility of the measurements had to 
be better than 0.5%. Therefore, care had to be taken 
to remove any factor that affects the data quality. The 
following somewhat unconventional procedure was per- 
formed to obtain reliable data. First, the detector slit 
was set wide enough (A26* ~ 4 degrees) to accept all the 
diffracted x-ray within a certain Bragg peak (Fig. 1). 
This is to minimize the detector non-uniformity effect as 
use of angle limiting devices such as Soller slits resulted 
in very unreliable data. However, some degree of non- 
uniformity still existed and it was necessary to make it 
sure that the same part of the detector is used for each 
scan. To ensure this, we scanned 29 at each temperature 
to locate the detector so that the centroid of the diffrac- 
tion peak is at the center of the detector. In this way, 
the detector is essentially tracking the temperature de- 
pendent motion (in 29) of the diffraction peak. As the 
last step, the 9 scans (rocking curves) were taken to mea- 
sure the diffraction peak intensity. The above procedure 
produced the most reliable data. 

Fig. 2 shows 9 scans as well as temperature depen- 
dent intensity plots measured as described above. Panel 
(a) shows normalized 9 scans of the (0,0,13) peak from 
YBCO. The peak position in 9 decreases as temperature 
increases, implying the increase of the c-axis lattice con- 
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FIG. 2: (a)Series of 9 scans of (0,0,13) peak from YBCO at 
different temperatures. The changes in total intensity as well 
as the peak position are seen. (b)Temperature vs. integrated 
intensity plot of the data. The data were taken with the 
temperature cycled and normalized at T=0. Also shown is 
the data for (0,0,11) peak. (c)Scans of (0,0,14) from Hgl212 
samples. (d)Temperature vs. intensity plot for (0,0,14) and 
(0,0,10) peaks. 



stant as expected. In addition to the decrease in the 9 
position, we note that the Bragg peak intensity decreases 
as the temperature increases (increased constant back- 
ground due to incoherent scattering was also observed 
but subtracted in the plot and analysis). No appreciable 
change in the diffraction line shape, that is, no q depen- 
dence is observed and thus the role of thermal diffuse 
scattering is not considered in the following discussions. 
Since the detector slit was wide open (in 29) to accept 
all the x-rays of the peak, each point on the 9 scan repre- 
sents Bragg peak intensity i.e., integration of a 29 scan at 
the given 9 value. Therefore, the integration of the area 
represents the sum of the intensities from all the grains 
of the sample. It is apparent from the figures that the 
peak intensity decreases as the temperature increases. 

To quantify the temperature dependent peak inten- 
sities, the peaks shown in Fig. 2a are integrated over 
the whole angular range. The results are normalized 
to a common linear extrapolated T=0 value and plot- 
ted against the temperature in panel (b). For monatomic 
systems, the intensity of a (0,0,Z) Bragg peak is expressed 
as I = I exp(— ^(z 2 )l 2 ) where the exponent is the DWF 
and (z 2 ) is the thermally averaged mean square mo- 
tion of the atoms in z direction—. For a classical Har- 
monic oscillator, (z 2 ) oc T and the intensity thus be- 
comes I = Ioexp(—Xl 2 T) oc 7 (1 — \l 2 T) where A is a 



constant. Therefore, one would expect linear tempera- 
ture dependence of the DWFs for T< l/\l 2 . The in- 
tensities in the figure show linear temperature depen- 
dence and no anomaly is found within the experimen- 
tal error—. Also shown in the panel is the temperature 
dependent intensity of the (0,0,11) peak from YBCO. 
The ratio of the DWFs of the (0,0,13) and (0,0,11) peaks 
is 1.28 ± 0.05 which is somewhat close to the expected 
value of 13 2 /11 2 = 1.4 for monatomic systems with pure 
thermal phonon effects. The measured temperature de- 
pendent intensities suggest normal behavior of the c-axis 
DWF without any anomaly. Panels (c) and (d) show 
data from Hgl212 samples. Other than the larger rock- 
ing width, the behavior is more or less the same with 
YBCO case, showing no sign of anomaly at T c within 
the experimental error. The deviation of (0,0,10) data 
from the linear behavior is related to the fact that, un- 
like other peaks, it was measured along with substrate 
peaks hence producing larger experimental errors, and is 
assumed to be extrinsic effect. 

Interpretation of the above results, however, is not as 
easy as the case for monatomic systems such as Si. The 
expression of a Bragg peak intensity contains in general 
exponential functions of the displacement parameters of 
the atoms in the unit cell. For a monatomic case, it 
could be reduced to a single exponential function of the 
displacement parameter (assuming the displacement pa- 
rameters for the atoms in the unit cell are the same) 
and the DWF can be interpreted with a relative ease. 
The form factors for polyatomic compounds however are 
different and the Bragg peak intensity is expressed by a 
combination of exponential functions. It is therefore very 
hard to extract the information on the dynamic displace- 
ments from DWFs. The other aspect is that the observed 
temperature dependence does not solely come from the 
dynamic displacements. Not only the phonons but also 
the changes in the static atomic positions within the unit 
cell affect the intensities. 

In spite of the difficulties, we can still extract useful in- 
formation from the results. Even though the Bragg peak 
intensity expression contains multiple exponential terms, 
it shows linear temperature dependence if the displace- 
ments are purely due to thermal vibrations. Therefore, 
the fact that experimental data show linear temperature 
dependence strongly suggests that the major contribu- 
tor to the DWF is the dynamic displacements as there 
is no intrinsic reason for such temperature dependence 
from the static displacements. Indeed, detailed neutron 
diffraction experiments on YBCO show that the mean 
square displacements are mostly dynamic 19 . We also 
note that heavy elements contribute more to the diffrac- 
tion intensities in our x-ray measurements due to their 
greater high Z sensitivity. To see the different contribu- 
tions, we list in Table 1 calculated fractions of the pla- 
nar oxygen and Cu-0 planar contributions to the total 
diffraction peak intensities of YBCO and Hgl212 using 
the published structural information— Note that the 
fractions do not add up to make the total due to the in- 
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Sample 


Fl (total) 


F c 2 (Cu-0 plane) 


Fl (planar O) 


YBCO 
Hgl212 


50.0 
1386 


2.58 
1.90 


1.00 
1.00 



TABLE I: Calculated F c 2 's for YBCO (0,0,13)and Hgl212 
(0,0,14) normalized to plane O values. Cu-0 plane and planar 
oxygen contributions have been calculated with the occupan- 
cies of the other atoms set to zero. 



terference but roughly show the insignificant Cu-0 plane 
contribution to the total intensity. It is thus conceivable 
that the phonon anomalies in the Cu-0 plane measured 
by EXAFS— may be buried under the contributions from 
the heavy elements. Based on these, it is reasonable to 
conclude that at least the heavy elements do not show 
anomaly at T c in their c-axis dynamic properties. 

On the more fundamental side, the anomalies in the c- 
axis lattice constants are generally relatively small com- 
pared to those along a- or b-axe a 4 i 6 i 7 . Therefore, it is 
possible that the c-axis phonon population anomaly may 
well be very small and undetectable. Looking at a- and fa- 
axis Dcbye- Waller factors is therefore essential. This may 
explain the absence of the anomaly in the Hgl212 data 
while measurement by EXAFS on powder samples of sim- 
ilar compound HgBa2Cu04 shows anomaly in the mean 
square dispacementsifi. Use of Rietveld refinement with 
powder diffraction at a brighter x-ray source could give 
an answer to this. Better yet, considering all the factors 
discussed above, similar experiments by neutron diffrac- 
tion may resolve most of the problems; the typical grain 
sizes of HTSC crystals and sample motion upon temper- 



ature change are insignificant compared to the neutron 
beamsize. This will allow us to investigate a- and b-axis 
DWFs by a diffraction technique. In addition, its sensi- 
tivity to low Z atoms will allow more accurate estimation 
of the oxygen contribution to the DWFs. There indeed is 
a strong indication that temperature dependent (1,2,12) 
neutron Bragg peak intensities of YBCO crystals have 
anomalies at the T,jiS. 

In conclusion, temperature dependent c-axis DWFs of 
YBCO and Hgl212 thin films are obtained by measur- 
ing the temperature dependent (0, 0, 1) x-ray Bragg peak 
intensities. The measured DWFs show linear tempera- 
ture dependence within the experimental errors, show- 
ing no anomalous effects near T c . Considering the fact 
that the x-ray diffraction intensity is much more sensitive 
to the heavy elements, the results suggest that there is 
no anomalous lattice dynamics of heavy elements along 
the c-axis. Further studies of a- and b-axis temperature 
dependent DWFs by neutron diffraction in combination 
with temperature dependent structural studies by EX- 
AFS may shed more light on the temperature dependent 
lattice dynamics, i.e., role of the phonons in HTSCs. 

We would like to thank F. Bridges and J. Arthur for 
helpful discussions, and H. Shin for statistical analy- 
sis. SSRL is operated by the DOE office of Basic En- 
ergy Research, Division of Chemical Sciences. The of- 
fice's division of Material Science provided support for 
this research. This work is supported (in part) by the 
Korean Science and Engineering Foundation (KOSEF) 
through Center for Strongly Correlated Materials Re- 
search (CSCMR) at Seoul National University. 



1 A. Nazarenko, and E. Dagotto, Phys. Rev. B 53, 2987 
(1996). 

2 T. Sakai, D. Poilblanc, and D. J. Scalapino, Phys. Rev. B 
55, 8445 (1997). 

3 H. You, U. Welp, and Y. Fang, Phys. Rev. B 43, 3660 
(1991). 

4 T. Asahi, H. Suzuki, M. Nakamura, H. Takano, and J. 
Kobayashi, Phys. Rev. B 55 9125 (1997). 

5 K. Kadowaki, F. E. Kayzel, and J. J. Franse, Physica C 
153-155, 1028 (1988). 

6 M. Braden, O. Hoffels, W. Schnelle, B. Biichner, G. Heger, 

B. Hennion, I. Tanaka, and H. Kojima, Phys. Rev. B 47, 
12288 (1993). 

7 V. Pasler, P. Schweiss, C. Meingast, B. Obst, H. Wiihl, A. 
I. Rykov, and S. Tajima, Phys. Rev. Lett. 81, 1094 (1998). 

8 C. Y. Yang, S. M. Heald, J. M. Tranquada, A. R. Mood- 
enbaugh, and Y. Xu, Phys. Rev. B, 38, 6568 (1988). 

9 C. H. Booth, F. Bridges, J. B. Boyce, T. Claeson, B. M. 
Lairson, R. Liang, and D. A. Bonn, Phys. Rev. B 54, 9542 
(1996). 

A. Lanzara, N.L. Saini, A. Bianconi, F. Due, and P. Bor- 
det, Phys. Rev. B, 59, 3851 (1999). 

C. Kittel, Introduction to Solid State Physics, John Wiley 
& Sons, Inc., New York, 1986. 



B. F. Cole, G.-C. Liang, N. Newman, K. Char, G. Za- 
harchuk, and J. S. Martens, Appl. Phys. Lett. 61 1727 
(1992). 

J. Z. Wu, S.L. Yan, and Y.Y. Xie, Appl. Phys. Lett. 74 
1469 (1999). 

a- or b-axis YBCO thin films do exist. However, these films 
are grown on LaSrGa04 substrates which have a cubic 
structure. The result is that the substarte peaks, which 
are much stronger, overlap with the peaks from the film, 
making those film sample useless for our experiments. 
This formula holds true only for monatomic systems. For 
polyatomic compounds, it is more complicated. See the 
discussion that follows later in the text. 
The fluctuations in Fig 2b and 2d are not from the sta- 
tistical shot noise in the counts but are mostly due to the 
extrinsic effects such as the thermal motion of the sample 
as discussed in the text. The size of the error bar based on 
our experience is 0.5% or less. All of the root mean square 
errors of the data are less than 0.5%. 

J. D. Jorgensen, B. W. Veal, A. P. Paulikas, L. J. Nowicki, 
G. W. Crabtree, H. Claus, and W. K. Kwok , Phys. Rev. 
B, 41, 1863 (1990). 

For Hgl212, published atomic positions of TlBa2CaCu207 
which is isostructural with Hgl212 were used. For 



TlBa 2 CaCu 2 07 structure, see B. Morosin, D. S. Ginley, P. 19 P. Schweiss, W. Reichardt, M. Braden, G. Collin, G. Heg 
F. Hlava, M. J. Carr, R. J. Baughman, J. E. Schirber, E. H. Claus, and A. Erb, Phys. Rev. B, 49, 1387 (1994). 

L. Venturini, and J. F. Kwak, Physica C 152, 413 (1998). 



